Prom local to critical fluctuations in lattice models: a non-perturbative 

renormalizat ion- group approach 
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We propose a new implementation of the non-perturbative renormalization-group (NPRG) which 
applies to lattice models. Contrary to the usual NPRG approach where the initial condition of 
the RG flow is the mean-field solution, the lattice NPRG uses the (local) limit of decoupled sites 
as the (initial) reference system. In the long-distance limit, it is equivalent to the usual NPRG 
formulation and therefore yields identical results for the critical properties. We discuss both a lattice 
field theory defined on a d-dimensional hypercubic lattice and classical spin models. The simplest 
approximation, the local potential approximation, is sufficient to obtain the critical temperature 
and the magnetization of the 3D Ising, XY and Heisenberg models to an accuracy of the order of 
one percent. We show how the local potential approximation can be improved to include a non-zero 
anomalous dimension 77 and discuss the Berezinskii-Kosterlitz-Thouless transition of the 2D XY 
model on a square lattice. 

PACS numbers: 05.70.Fh,05.10.Cc,05.70. Jk 



I. INTRODUCTION 

Many models of statistical physics and condensed mat- 
ter are defined on a lattice. The phase diagram of a lat- 
tice model usually depends on the lattice type, the range 
and strength of interactions, as well as other details of the 
Hamiltonian. On the other hand, the precise knowledge 
of the Hamiltonian is often not necessary to understand 
the long-distance behavior of the system, in particular 
the universal critical properties near a second-order phase 
transition. In this paper, we describe an approach based 
on the non-perturbative renormalization group (NPRG) 
which captures both local and critical fluctuations in lat- 
tice models and therefore describes universal and non- 
universal properties. In particular, we can calculate the 
critical exponents, the transition temperature and the 
magnetization in classical spin models. 

The NPRG approach has been successfully applied to 
many areas of physics and in particular to the study of 
critical phenomena [l}-^. It has recently been extended 
to lattice models [4]. The strategy of the NPRG is to 
build a family of models indexed by a momentum scale 
parameter fc, such that fluctuations are smoothly taken 
into account as k is lowered from a microscopic scale A 
down to 0. In practice, this is achieved by adding to 
the Hamiltonian (or the action) a "regulator" term AHk , 
which vanishes for fc = 0, and computing the correspond- 
ing Gibbs free energy (or effective action in the field the- 
ory terminology) F^. The initial value Ai/^ is chosen 
such that in the reference system defined by the Hamil- 
tonian H -f Ai/A all fluctuations are effectively frozen. 
The determination of Fa then reduces to a saddle-point 
(mean-field) calculation. The Gibbs free energy Ffc=o we 
are eventually interested in is obtained from that of the 
reference system by solving a RG flow equation. The 
latter cannot in general be solved exactly (even numeri- 
cally) and one has to resort to some approximations. The 



approximate flow equation must be sufficiently accurate 
(and yet tractable) to provide a good approximation of 
the state of the system. 

In some cases however, the mean-field solution is too 
far away from the actual state of the system to provide a 
reliable initial condition for the NPRG procedure. An ex- 
ample is provided by the localization transition between a 
Mott insulator and a superfluid in lattice boson systems. 
The two-pole structure of the local (on-site) propagator 
is crucial for the very existence of the transition. This 
structure is however impossible to reproduce using a RG 
approach starting from the mean-field (Bogoliubov) ap- 
proximation. This prevents a straightforward generaliza- 
tion of recent NPRG studies of interacting bosons 
to lattice models such as the Bose-Hubbard model. 

We therefore propose a new NPRG scheme for lattice 
models where the reference system corresponds to the (lo- 
cal) limit of decoupled sites. As an expansion about the 
local limit, the lattice NPRG is reminiscent of Kadanoff's 
idea of block spins fT2], although the way intersite inter- 
actions are progressively introduced when lowering the 
momentum scale k makes it significantly different from 
a real-space RG. In the long-distance limit, the lattice 
NPRG is equivalent to the usual NPRG formulation and 
therefore yields identical results for the critical proper- 
ties. 

The possibility to start from a reference system which 
already includes short-range fluctuations has been recog- 
nized before, and was used by Parola and Reatto in the 
Hierarchical Reference Theory of fluids [l^ , an approach 
which bears many similarities with the lattice NPRG. An 
important aspect of the lattice NPRG is that it is formu- 
lated in the field theory language commonly used in the 
NPRG approach. Its relation to the standard NPRG 
formulation is therefore obvious, and many of the ap- 
proximate solutions of the flow equation satisfied by Ffe 
proposed previously also apply to the lattice case. 
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In Sec. |TT]we introduce the lattice NPRG for a lattice 
field theory. We first recall the "standard" NPRG ap- 
proach to lattice models Q , and then show that the lat- 
tice NPRG scheme merely results from a different initial 
condition while the long-distance (small k) behavior of 
the effective action remains the same. As an applica- 
tion, we derive a lattice field theory from the Ising model 
and compute the transition temperature in the local po- 
tential approximation (LPA). In Sec. IIIIl we apply the 
lattice NPRG to classical spin models without deriving 
first a lattice field theory. We find that the LPA is suffi- 
cient to obtain the critical temperature and the magneti- 
zation of the 3D Ising, XY and Heisenberg models to an 
accuracy of the order of 1 percent. We also discuss an im- 
provement of the LPA (known as the LPA') which yields 
a non-zero anomalous dimension 77. In Sec. IIVI we use 
the lattice NPRG to calculate the Berezinskii-Kosterlitz- 
Thouless (BKT) [ij, [T^ transition temperature of the 
2D XY model on a square lattice. 

II. LATTICE NON-PERTURBATIVE RG 

We consider a lattice field theory defined on a d- 
dimensional hypercubic lattice, 

q r 

where {r} denotes the TV sites of the lattice. For sim- 
plicity, we consider a one-component real field ipr- y^q — 
(^r is the Fourier transformed field. The 
momentum q is restricted to the first Brillouin zone 
] — IT, tt]'^ of the reciprocal lattice. In the thermodynamic 
limit (A^ — > 00), 




The potential [/q is defined such that eQ{q = 0) ~ 
but is otherwise arbitrary. eo(q) — eoQ^ for q — and 
maxqeo(q) = e™'^^- The lattice spacing is taken as the 
unit length. 

To implement the RG procedure, we add to the Hamil- 
tonian ([T]) the regulator term 

^Hk[^] = ^Y.^-^^k{ci)iPq. (3) 
q 

Throughout the paper, we take 

i?fc(q) = (efe-eo(q))0(efe-eo(q)) (4) 

(f'k = eo^^), which is adapted from Ref. to the lat- 
tice case. The cutoff function i?fc(q) leaves the high- 
momentum modes (eo(ci) > unaffected and gives a 
mass to the low-energy ones (their effective (bare) dis- 
persion satisfies eo(q) -I- i?fe(q) = e^). 



In the presence of an external field, the partition func- 
tion reads 

Zk[h] = J 2?[(^]e-f^M"^ff''M+5:.''r'^r (5) 

and the order parameter is given by 

din Zk[h] 

</)r = {(fir) = — . (6) 

The so-called average effective action, 

Tk[cb]^ -lnZk[h]+Y,hrcbr- AHk[cb], (7) 

r 

is defined as a modified Legendre transform which in- 
cludes the explicit subtraction of AHkicf] [2]. It satisfies 
the exact flow equation [l| 

q 

as the energy scale efe is varied. Fj, [0] is the second- 
order functional derivative of F^ [</<]. Since Rk=o{(l) ~ 0, 
Ffe=o [0] coincides with the effective action of the original 
model (H]). 

A. Standard NPRG scheme 

In the "standard" NPRG approach to lattice models Q , 
the initial value A of the momentum scale k is chosen such 
that eA is much larger than all characteristic energy scales 
of the problem. In this limit, all fluctuations are frozen 
and mean-fleld theory becomes exact: T\[(j>] = H[(j)]. The 
first part of the RG procedure, when varies between eA 
and ek^^ — e™^^ is purely local, since the effective (bare) 
dispersion eo (q) + Rk (q) — Cfc remains dispersionless for 
all modes. Only for k < fcjn does the intersite coupling 
start to play a role. When fc ^ 1, i.e. when I/fc is 
much larger than the lattice spacing, the lattice does not 
matter any more. This result is a direct consequence 
of the structure of the flow equation; the dkRk term in 
^ implies that only modes with |q| ^ A; contributes 
to 9feFfc. When fc ^ I, one can therefore approximate 
eo(q) — eoQ^i and one recovers the flow equation of the 
continuum model obtained from ([T]) and ([3]) by replacing 
eo(q) by eoq^ 0. 

B. Lattice NPRG scheme 

In the lattice NPRG, we start the RG procedure from 
fc = fcin, i.e. we bypass the initial stage of the flow fcin < 
fc < A where the fluctuations are purely local (Fig. [T]) . 
The average effective action Ffcjjj[(/)] is no longer given by 
the microscopic Hamiltonian H[(l)\ (since the mean-field 
solution is not exact for the Hamiltonian H + AHk^^ ) but 
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Sec. Ill C[ this is the case for classical models even within 
simple approximations. In practice however, one often re- 
lies on an approximate solution of the flow equation. The 
NPRG lattice scheme is preferable whenever the (approx- 
imate) flow equation gives a poor description of Vk^^ [0] 
starting from the mean- field result rA[0] = As dis- 

cussed in the Introduction, this is to be expected in quan- 
tum models (such as the (Bose-)IIubbard model) where 
on-site (quantum) fluctuations make the local limit non- 
trivial. Finally we point out another advantage of the 
lattice NPRG; it enables to study classical spin models 
without first deriving a field theory (see Sec. IIII[) . 



FIG. 1: (Color online) Initial effective dispersion eo{q) + Rk {q) 
in the standard (left panel) and lattice (right panel) NPRG 
schemes (d = 1 and eo(f?) = 2eo(l — cos q)). The green dashed 
line shows the bare dispersion eo(<?). 

[k = fcui) > fc > 0) (fc = 0) 




FIG. 2: (Color online) Effective (bare) dispersion eo{q)~\-Rk[q) 
for k = fcin, fcin > > and fe = with the cutoff function 
(|4]) (d = 1 and eo[q) = 2eo(l — cosg)). The green dashed line 
shows the bare dispersion eo(?)- 



C. Application to the Ising model 

We consider the Ising model, 

H = - J/3 SrSr' (S'r = ±1) (10) 

defined on a d-dimensional hypercubic lattice (/3 = 1/T). 
(r, r') denotes nearest-neighbor sites. To apply the 
NPRG approach, one possibility is to first derive a field 
theory (another, more natural, approach is described in 
Sec, imp . To this end, one considers the Hamiltonian 



- ^ SrAl'^l Sr' , 



(11) 



its computation reduces to a single-site problem which 
can be easily solved numerically (and even analytically 
in some models). In principle, the NPRG scheme can be 
defined with any cutoff function provided that the latter 
satisfies the initial condition [13] 



i?fc.„(q) --eo(q)+C 



(9) 



ensuring that the sites are decoupled (the limit C — )■ cx) 
corresponding to the standard scheme) . The choice C — 
Efcin = made in ^ is however very natural, since 
it allows to set up the RG procedure for k < fcin in the 
usual way, i.e. by modifying the (bare) dispersion of the 
low-energy modes eo{q) < eu without affecting the high- 
energy modes (Fig. [2]). 

For fcin > fc > 0, the effective coupling in real space 
(defined as the Fourier transform of eo(q) + -Rfc(q)) is 
long-range and oscillating. The oscillating part comes 
from the behavior of -Rfc(q) for eo(q) ~ ffc- Although 
the lattice NPRG is based on an expansion about the 
local limit, it markedly differs from Kadanoff's real-space 
RG [iJI in the way degrees of freedom are progressively 
integrated out. 

The standard and lattice NPRG schemes thus differ 
only in the initial condition. Both schemes are equiv- 
alent for fc < fcin when the flow equation ([5]) is solved 
exactly for fc > fcin in the standard scheme. As shown in 



which differs from that of the Ising model only by the 
additive constant — /i/3iV. The matrix A'^^ is diagonal in 
Fourier space with eigenvalues 



Ap(q) = 2^f J^cosq^, 



(12) 



For II > Jd the matrix A^'^' is positive (A^(q) > Vq) 
and can be inverted. We can then rewrite the partition 
function of the Ising model using a Hubbard-Stratonovich 
transformation. 



{Sr} 



oc 



ip^f 111 cosh iff 



(13) 



We thus obtain a lattice field theory with the Hamilto- 
nian 



H, 



q 



AM(q) A^(0) 



2A^(0) 



In cosh (^r 



(14) 
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FIG. 3: (Color online) Left panel: Function h{(j>) obtained 
from the numerical solution of (|17|) . Right panel: Effective 
potential Uki^{(j>) [Eq. (HD)]. = 5, T = 4.48 J and d = 3. 



Rescaling the field, we can cast the Hamiltonian in the 
form ([T]) with 



eo(q) = 2d{Jd + ^i) 



l-7q ^ 



(15) 



max 




(7q — d^^^^cosq^). The (bare) dispersion eo(q) in- 
cludes long-range interactions, and eo(7r, tt, • • • ) — e, 
diverges for fi Jd. In the limit — > oo long-range 
interactions are suppressed and eo(q) ~> 2(i(l — 7q). 

We are now in a position to apply the NPRG ap- 
proach. In the standard scheme (Sec. Ill A[) . the initial 
value Tj^[(f)] = -ff [0] of the effective action is defined by 
the Hamiltonian ([1]), with eo(q) and Uq given by 
In the lattice NPRG scheme, the initial value Uki„ of the 
effective potential has to be computed numerically. One 
has Zki^[h] = Or ^fcin(^r), where 



(16) 



is the partition function of a single site in an external 
field h. The relation between (f>r and is obtained from 
the equation 



d 
dhr 



In Zk.Jhr), 



(17) 



which has to be computed and inverted numerically 
(Fig. [3|). The initial value of the average effective action 
then takes the form 



r r 
r q 

where 



I'qj 



(18) 



= -lnZfc,„(/i)+ 



(19) 




0.02 0.04 



0.06 



Inc 



FIG. 4: (Color online) Left panel: Derivative U'i^[p) of the 
effective potential for various values of k ranging from A to 
fcin. The initial value Ua{p) ~ Uo{p) is given by ([TSj. The 
red points show the solution 17^.^ (p) directly obtained from a 
numerical solution of the single-site partition function Zk^,^ (h) 
[Eq. (|16|) ]. Right panel: Relative error e vs a = (A/fcin)^. 
= 5, T = 4.48J and d = 3. 



is the effective potential (Fig. [3]) . 

In the local potential approximation (LPA), one ne- 
glects the /c-dependence of the dispersion so that 

TfcM = J2 + '/'-qeo(q)0q. (20) 



Here and in the following, we consider [/fe as a function 
of pr = </'?/2. From ([5]), one deduces 



dkUkip) = - 



afei?fe(q) 



2 eo(q) + i?fc(q) + Ulip) + 2pU'^{p) ' 



(21) 



With the cutoff function (U), equation (PT|) reduces to 



kdkUkip) 



6,. + C/^(p) + 2pC/^'(p) Jq 
The integral over q can be rewritten as 

/ 9{ek-eoiq)) = f ' deV{e) 

Jq Jo 

where [Til 



e{ek - eo(q))- 



P(6) 



£o(q)). 



(22) 
(23) 

(24) 



In the local fiuctuation regime k^^ < k < A, H[(p] + 
AHk [(fil is a local Hamiltonian (no intersite coupling) . It 
follows that both — In Zk [h] and its Legendre transform 
reduce to a sum of single-site contributions. The LPA 
for the average effective action Tk[(l)] is therefore exact. 
We have computed the derivative C/^ of the effective po- 
tential for various values of a = ^A/f-o^^^ = (A/kin)^. (In 
practice, it is easier to solve for U'l. than Uk-) The results 
are shown in Fig. [5] for a = 1000 (we comment on the 
numerical method in Sec. IIII A[) . As k decreases from 
A to kin, the potential U'i.{p) converges towards the ex- 
act solution obtained from the numerical computation of 
the local partition function Zki^{h) [Eqs. (|16I19I) ]. The 
relative error is shown in the right panel of Fig. [H It 
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decreases as 1/a, in agreement with the fact that the va- where 
hdity of the mean-field (saddle-point) approximation to 
Z\[h] for large A is controlled by R]^^ ^ e^^ ^ a^^. We 
therefore conclude that the standard and lattice NPRG 
schemes are equivalent in the LPA for classical lattice 
field theories. 

The critical temperature is obtained from the diver- 
gence of the susceptibility x = l/t^fe=o(/' ~ ^) (°^' 
equivalent!;^ the divergence of the correlation length 
^ = ^) [ig|. For a = 1000 and d = 3 one finds that 
Tc ~ 0.747 is independent of fi and in very good 
agreement with the "exact" result T^'''"'^ ~ 0.752 T^'^p 
obtained from Monte Carlo simulations l20i. 



^hS 



2 cosh(/i) 



(30) 



is the partition function of a single site in an external 
field h. The magnetization at site r, 



d 

rrir = — — \nz(hy.) = tanh(ft,r), 
ohr 



(31) 



varies between -1 and 1. Up to an additive constant, we 
obtain 



(32) 



III. CLASSICAL SPIN MODELS 

In this section, we show that the lattice NPRG can be 
applied to classical spin models without first deriving a 
field theory. For simplicity, we consider the Ising model 
on a c?-dimensional hypercubic lattice [Eq. ([TU])]. In the 
presence of an external field and a regulator term AHk, 
the partition function reads 

{Sr} 

^ ^ g-^EqS'_q[eo(q)-2€od-H«fc(q)]Sq+j:,'i,S,^ 



{Sr} 



where 



eo(q) 



d 

2eoE(l - cosq^.). 



(25) 



(26) 



with eo = J/T. Since 5^ = 1, the term 2deo in ([25]) 
contributes a constant term to the Hamiltonian and can 
be omitted. The magnetization at site r is given by 



d\nZk[h] 
dhr 



and the average efl:ective action is defined by 



Ffc [m] = - In Zk [h] 



AHklm] 



(27) 



(28) 



The standard NPRG scheme cannot be used, since the 
partition function is not expressed as a functional in- 
tegral over a continuous variable. A "regulator" term 
efe S-qSq = Efc J2r "^r = would Only add a con- 
stant term to the Hamiltonian. On the contrary, there 
is no difficulty to apply the lattice NPRG scheme. With 
the cutoff function (U), one has 



and the effective potential [2l| 



C/fc.„(p) = ^ln(l-2p) 



2/9 atanh( V 2p) - Meop, (33) 



where = mj/2. Uk,^{l/2) = ln(2) - 2c?eo is finite 

but J/^Jp) iln(l ~2p) diverges for p -> 1/2. This 

divergence suppresses the propagator 



1 



eu+ RM) + U'k{p) + '^pU'ap) 



(34) 



appearing in the flow equation ((2T|) and therefore the 
fluctuations corresponding to a large magnetization. 

A comment is in order here. We have followed the usual 
convention to deflne the average effective action F^ as a 
modifled Legendre transform which includes the explicit 
subtraction of AiJfe[m] [Eq. O] [2]. The deflnition of the 
average effective action F^ is of course arbitrary provided 
that Ffc=o corresponds to the true Legendre transform of 
the original model. Since 



UuM = p{l-Ade^) + 0{p') 



(35) 



for p — )• 0, we find that the initial transition temperature 
is determined by 1 = 4deo, i.e. Tc^'"'^ — MJ, which differs 
from the mean-field transition temperature = 2dJ 
by a factor of 2. Ffej^[m] assumes a mean-field treat- 
ment of the intersite coupling term i 5_qeo(q)S'q as 
in the usual mean-field approach to the Ising model. 
The discrepancy between Ti*"'"^ and T^^ comes from 
the fact that 5 X]q '^-q^o(q)5'q includes a local term 
■i 2deo>5'r = NdeQ. The latter contributes a mere con- 
stant to the Hamiltonian, but is considered at the mean- 
field level in the average effective action where it gives a 
term — deg rn^. To make contact with the usual mean- 
field theory, we consider the effective average action 



Tk [m] = Tk [m] 



r, r 



(36) 



Zkjh] 



{S.} 



'2deoN 



(29) 



and the corresponding effective potential 
Uk{p) = Uk{p) + pRk{r,r). 



(37) 
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rpMF 

^ c 




rriNPRG 

c 


Ising 3D 


6 


4.51 


4.48 


XY 3D 


3 


2.20 


2.18 


Heisenberg 3D 


2 


1.44 


1.42 



TABLE I: Critical temperature T^^^'^ obtained in the LPA 
compared to the mean-field estimate and the Monte- 

Carlo result T^""""* All temperatures are in unit 

of J. 




FIG. 5: (Color online) Transition temperature ri*^' obtained 
from the effective potential Uk (green points) and Uk (red 
points). 



rfc[r7i] differs from the true Legendre transform only by 
non-local terms. The initial value 



(38) 



reproduces the mean- field result Tc'^"'^ — . Again we 
stress that F^. and lead to the same fc = results and 
in particular to the same critical temperature. 



A. Local Potential Approximation 

We have solved the equation for [/(.(p) (see Eq. ([H])) 
numerically using Euler's method with a typical RG time 
step M = -10-4 {t = ln(fc/fci„)). The function f7((p) is 
discretized with a few hundreds points in p. The con- 
vexity of the potential (see Fig. [H]) makes the numerical 
resolution difficult below Tc (in particular at low temper- 
atures) and there is a tendency to numerical instability 
for large values of \t\. However, the value pQ for which 
U'f.{p) = usually converges before instability problems 
arise. 

The critical temperature is obtained from the criterion 
U^^^ip = 0) = (Sec. EC]). One finds ~ 0.747 T^^^^ 
for the three-dimensional Ising model, in very good agree- 
ment with the results of Sec. Ill CI and the "exact" result 
ycxact _ 0.752 T^'^^ obtained from Monte-Carlo simula- 
tions ■ We have obtained a similar accuracy for the 
critical temperature of the XY and Heisenberg models in 
d = 3 (Table H. 




FIG. 6: (Color online) Left panel: Potentials [7^=0 (p) (red 
solid line) and C/^. (p) (green dashed line) at criticality (T = 
Tc) in the LPA id = 3). Right panel: Potential U'k^o{p) for 
T = Tc (red dotted line), T = l.OSTc (green dashed line) and 
T = 0.95 Tc (blue solid line). 



Fig. [5] shows the transition temperature Tc^'^ of the 
three-dimensional Ising model deduced from the effective 
potentials Uk and Uk [Eq. ([57|) ]. As k decreases, Tc^^ con- 
verges rapidly towards the actual transition temperature 
Tc — Tc^~'^\ This result is due to the fact that all de- 
grees of freedom contribute more or less equally to the 
thermodynamics. Once k <^ fcin, thermodynamic quan- 
tities are therefore obtained with a reasonable accuracy. 
This also explains why the LPA, which does not correctly 
describe the long-distance limit of the propagator when 
T ~ Tf,, is remarkably successful in computing the tran- 
sition temperature and other thermodynamic quantities. 

Fig.[n]shows the derivative U'f.{p) of the effective poten- 
tial for k = /cin and fc = at criticality (T = 0.747 T^'^^), 
as well as J7^=o(p) ior T = T > Tc and T < Tc. 
In the latter case, we find UI^^q{p) = and therefore 
Uk=o{p) — const for p < po, where po = ?Tig/2 deter- 
mines the actual magnetization mg of the system. This 
result is a consequence of the convexity of the potential 
in the low-temperature phase, a property which is known 
to be satisfied in the LPA 

In Fig. [71 we show the uniform susceptibility \ = 
l/Ul^^f^{p = 0) in the high-temperature phase [l^, as 
well as the magnetization below Tc with the Essam- 
Fisher approximant [l^l- We find the critical exponents 
1/ = 2/3 = 7/2 ~ 0.64 - 0.65 (with r? = in the LPA), 
in agreement with the known result in the LPA with the 
cutoff function © [H, [111 . 

Not surprisingly, the LPA is not as accurate in two 
dimensions. For the 2D Ising model, we find Tc ~ 
0.48 T^^, to be compared with the exact resuh T^^'^'^*- = 
2J/ln(l + \/2) ~ 0.567 Tc'^'" l2i|. 



B. Renormalization of the spectrum 

A natural generalization of the LPA includes a renor- 
malization of the amplitude of the spectrum. We there- 
fore consider the Ansatz 



TfcH =^[/fc(pr) + i^Afceo(q)T 



_q//i,q. 



(39) 
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FIG. 7: (Color online) Top panel: Uniform susceptibility x 
(red points) in the high-temperature phase and a fit x oc 
(T — Tc)~^ with 7 = 2j/ ~ 1.30. Bottom panel: Magnetization 
m in the low-temperature phase (red points) and the Essam- 
Fisher approximant [24'| (solid line). The inset shows a fit to 
m « (Tc - Tf with P = v/2~ 0.32. 



This approximation can be seen as the first step of a cir- 
cular harmonic expansion of the renormalized dispersion 
e(q) t^]. Since 

ri')(q; p) - A,6o(q) + U^P) + 2pC/^'(p) (40) 

in a uniform field (p^ — %/2p, we can define the renormal- 
ized spectrum amplitude by 

1 



Co 



(41) 



•where r and r' are nearest neighbors. The amplitude 
Ak = Ak{po^k) should be understood as the first term in 
the expansion of the function 



Akip) = Ak{pQ,k) + A[}\po^k)ip - Po,k) 



(42) 



about the minimum po.fc of the effective potential Uk{p)- 
Another possible definition of the spectrum amplitude 
is Ak = Ak{pa,k), "where po,k is the minimum of Uk{p) 
[Eq. ([57| ]. The flow equation for Ak follows from ([5]) and 

dSD, 



dkAk = ^ 
eo 



1 - 



1 - 



eo(q) 

£o(p) ' 
2den 



dkRk{ci)G{cif 
G(p), (43) 



where 73 = y^2p^kW^{pa,k) + 'ipa,kU',l' {po,k)]- The flow 
equation for Uk{p) is identical to (PT|) with eo(q) replaced 
by Afceo(q)- 




FIG. 8: (Color online) Left panel: Afc(po,fe) (red solid line) and 
Ak{po,k) (green dashed line) vs ~t = ln(fein/fc). Right panel: 
-fcflfc In Afe(po,fc) (red solid line) and ~kdk\riAk{po,k) (green 
dashed line) vs —t. All curves are obtained at criticality (T = 
Tc) and for d = 3. 



Numerical results for Ak{po,k) and Ak{po,k) are shown 
in Fig. [51 In both cases, Ak varies when fc ~ 1, in agree- 
ment with the expectation that the amplitude of the har- 
monic coii{nqu) should vary when fc ^ 1/n [4]. The vari- 
ation of Ak is moderate (from 1 to 1.27 for Ak{po,k) and 
from 1 to 1.19 for Ak{po,k)) and weakly affects the crit- 
ical temperature which remains within a few percents 
of the exact result: = 0.74 T^^'^ with Ak{po,k) and 
Tc ^ 0.72 T^^P with Ak{po,k)- We expect that the inclu- 
sion of additional higher-order harmonics in the spectrum 
would give a better estimate of Tc- 



C. LPA' 

When fc ^ 1, we can approximate eo(q) by eoQ^- In 
this regime, a simple improvement over the LPA (known 
as the LPA') consists in including a field renormalization 
factor Zk so that the renormalized dispersion e(q) is given 
by Zfceo(q) — ZkCoq^. The LPA' can be generalized to all 
values of k by writing the average effective action as [4] 



PfcH = ^ L/fc(/?r) + I ^m_qZfceo(q)mq. 



(44) 



Although this Ansatz is formally similar to ((39l) . Zk 
should not be confused with the amplitude Ak introduced 
in the preceding section. Zk is computed from the 0{q^) 
part of the spectrum, 



1 d 



(45) 



1 d 
— hm 7— . 

and therefore receives contributions from all harmonics. 
The LPA' can be justified when fc ~ 1 by noting that 
in this limit the renormalization of the spectrum is weak 
(Zk ~ 1), so that the approximation e(q) ~ Zfceo(q), 
valid for small q, is expected to remain approximately 
valid in the whole Brillouin zone Q. Nevertheless, be- 
cause short-range fluctuations are important for the ther- 
modynamics (Sec. lIII A]) , the LPA' might lead to a slight 
deterioration of the value of Tc obtained in the LPA. As 
in the preceding section, we can compute Zk either from 
the minimum of Ukip) (as in (j45p ) or from the minimum 
Po.fe of C/fe(p). 
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To obtain a fixed point when the system is critical, one 
should redefine the cutoff function, 

i?fe(q) = Zk{ek - eo(q))e(efc - eo(q)), (46) 

and introduce the dimensionless variables 

p - Zkk~^ekP. Ukip) = k-''Uk{p). (47) 

This change of variables cannot be done at the beginning 
of the flow if one works with a given range of p values. 
This would indeed correspond to a smaller and smaller 
range in p, whereas a good determination of the critical 
temperature requires to consider the window < p < 
1/2. To circumvent this difficulty, we define 

p = Zkk-''ekgi{k)p, Ukip) = k~''g2{k)Uk{p), (48) 
where 

' k^iZkCk)-^ if k » fc„ 



9i{k) 



92{k) = 



1 



if fc < fee, 



fc'' if fc > kc 
1 if fc < fee 



(49) 



For k ^ kc, p and Uk are equal to p and Uk, whereas 
they coincide with the dimensionless variables (j47|) when 
k <^ kc- The momentum scale kc will be determined 
below. In practice, we take 



gi{k) = [x + k-''Zkek{l-x)]-\ 
g^{k)^[x + k-''{l-x)r\ 



(50) 



where x — e^'^^^^"^ with n > d. The parameter n fixes 
the size of the crossover region k kc between dimen- 
sionful and dimensionless variables (Fig. [9]). 

Fig. [To] shows the RG time —t{p) beyond which the 
potential U'^^p) differs from [/^^q(p) by less than 10"'^ 
{d = 3). (The RG time t is defined by fc = fcine*.) 
This precision (10~^) is sufficient to determine the criti- 
cal temperature to an accuracy of the order of 1 percent. 
As expected, the potential converges to its asymptotic 
value faster for large values of p. In Fig.[Tni the red solid 
line shows the maximum value Pmax = time r if 

one works with a fixed window p € [0, 1/2] [Eq. (j47|) ] (we 
neglect Z^^c^ ~ 1 in the calculation of Pmax)- Thus we 
see that a natural choice for tc = ln(fcc/fcin) is tc — —2.5, 
since for times —t > —tc the window [0,/9max] becomes 
larger than the range of p values for which the potential 
U'f.{p) has not converged to its asymptotic value yet. In 
practice, we verify that our results are independent of the 
precise choice of kc and n. 

The change of variables (pS)) leads to the flow equations 



kdkUk 



-d+'-^Wk 



+ 



.92 

2 + % 

92 



kdkgi \ 



pu'k 



91 ) 

[{2^r^k)h+Vkl2] (51) 




1 2 3 4 5 6 7 
— t = ln(fcijj/fc) 

FIG. 9: (Color online) The function /2(fc) = k-'^g2{k) vs 
— t = — ln(fc/fcin) for n — 6 and tc = —3 {d — 2). The green 
dashed line corresponds to g2{k) = 1. 




0.2 0.3 
P 



0.4 



0.5 



FIG. 10: (Color online) RG time ~r{p) beyond which the 
potential U'j.{p) differs from UI.-q{p) by less then 10~^ (blue 
dotted line). The red solid line shows p — (see text). 




FIG. 11: (Color online) Left panel: Zk vs ~t = ln(fcin/fc) 
obtained with po.fc (red solid line) and po.fe (green dashed line) 
(d = 3). Right panel: Anomalous dimension 77^ = —kdkln Zk 
obtained with pa.k (red solid line) and po.fc (green dashed line). 



and 

^k^A^p[W'^ + 2~pU'lf 

92 



l3 + Sd.2/{8Tr) 



l + i^iU', + 2pU'^) 



(52) 



where 



eo(q), 



'(efe - eo(q)), 



u2-d 



£o(q)^g.£o(q) - [^g.£o(q)]^ 
eo(q)' 



(53) 



Kefc-eo(q)) 



(see Appendix [XI) . In Eq. ([5^ . the rhs should be eval- 
uated at po.fe or pa,k- For d — Z and n > 3, the nu- 
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merical solution of the flow equations (|51I52I) is stable 
for —tc > 2. The results described below are obtained 
for n = 4 and = -3. We find = O.ST^'^ when 
Zk is defined with respect to the minimum po.fe of Uk, 
and Tc — 0.74 T^'^^ if we use the minimum po.fe of Uk- 
The result is not as accurate as in the LPA (as an- 
ticipated above; see the discussion following Eq. (|45l) ). 
Nevertheless, with po.fe (the only case we discuss in the 
following), it remains within 2 percent of the exact re- 
suh T^''^'^' 0.752 T^'^^. The flow of Zk is shown in 
Fig- [HI For —t<l,Zk does not differ significantly from 
Ak (Fig. [5]). In this regime, only the first harmonic (i.e. 
cosg^) is expected to vary and therefore contribute to 
Zk- For — i > 1, the renormalization of higher-order har- 
monics makes Zk deviate from Ak- While Ak saturates 
to ~ 1.19, Zk ^ k^^ diverges with an exponent given 
by the anomalous dimension rj* = liiak^Q rjk- rj* ~ 0.1 
is a poor estimate of the exact result rj* ~ 0.036 but 
agrees with previous estimates based on the LPA' [4]. A 
P dependence of Zk is expected to improve the value of 



P{PQ.k) = kdkPo.k nearly vanishes for po,fc > p^i which 
reflects the existence of a line of quasi-fixed points. In 
this low-temperature phase, after a transient regime, the 
running of the stiffness po,fe becomes very slow, which 
implies a very large, although not strictly infinite, corre- 
lation length ^. The anomalous dimension rjk depends on 
the (slowly varying) stiffness and takes its largest value 
when the system crosses over to the disordered regime 
(po./c ~ Po ^'^"i ^ ^ ^~^)- When /5o,fe < Pq, the essen- 
tial scaling ^ ~ Q°-/{Po~Po.k) of the correlation length is 
reproduced [37] . 

In this section, we apply the lattice NPRG to the two- 
dimensional XY model defined by the Hamiltonian 



H 



Si' ■ Sr 



(54) 



where Sr = (cos6'r, sin^r) is a 2D classical spin of unit 
length. Up to a multiplicative constant, the initial value 
of the partition function ^fei„[h] = Hr^C'^r) (^r = |hr|) 
is determined by the partition function of a single site. 



D. Comparison with HRT 



zih) 



dO 
2^' 



loih). 



(55) 



Our approach bears similarities with the Hierarchical 
Reference Theory (HRT) of fiuids (H-IH (for a review, 
see Ref. [13 )• The HRT is based on an exact treatment of 
short-distance (hard-core) interactions supplemented by 
a RG analysis of long-range interactions. The HRT also 
applies to classical spin models: as in the lattice NPRG, 
it starts from the local theory (decoupled sites) and takes 
into account the intersite coupling in a RG approach. Al- 
though the final results are very similar to those we have 
obtained in Sees. IIII Al and IIII Bl the HRT nevertheless 
differs from the lattice NPRG in some technical aspects, 
e.g. the choice of the cutoff function and the way de- 
grees of freedom are progressively integrated out (Tsl . [33| . 
Because it was first developed in the context of liquid 
state theory, the connection between HRT and the more 
standard formulation of the NPRG |i2j] is not always ob- 
vious (for a discussion of the relation between HRT and 
RG, see Refs. [H,!!!!). By contrast, the lattice NPRG is 
formulated in the usual language of statistical field the- 
ory. The various improvements over the LPA known for 
continuum models can then be easily implemented in the 
lattice NPRG. In Sec. IIII CI we have discussed one of 
these improvements, the LPA', which allows to compute 
the anomalous dimension 77. 



IV. BKT TRANSITION IN THE 2D XY MODEL 

The NPRG approach to the continuum 0(2) (linear) 
model reproduces most of the universal properties of the 
BKT transition in two dimensions [36, 37]. In particular, 
one finds a value pg of the dimensionless order param- 
eter (the spin-wave "stiffness") such that the /3 function 



The magnetization points along the applied field with an 
amplitude 



m(/.).-lnz(;^) = ^. 



(56) 



where Io{h) and Ii{h) are modified Bessel functions. 
Contrary to the Ising model, the function h{m) obtained 
by inverting (j56p must be computed numerically. For 
k = kirn the average effective action takes the form 



Ffc^Jm] = ^ Uk,M + IY1 «o(q)m_q • niq (57) 



(pr — m^/2) with 



C^fc.n (P) = ~'inz{h) + h^2p ~ 8eop 



(58) 



Expanding Io{h) and Ii{h) for small h, one finds m = 
h/2 + 0{h^) and in turn 



C/fe,„(p) = 2p(l-4eo)-fO(p2). 



This yields the transition temperature Tc'^'"' = 4 J, which 
differs from the mean- field result = 2 J for reasons 
explained in Sec. IIIII 

In the LPA', the fiow equations read 



(59) 



kdkUk= {-2+'^^Uk+[Vk 

92 



^ ) pK 

92 



X [(2 - r]k)Ii + Vkh], 



(60) 
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FIG. 12: (Color online) Flow trajectories {po,ri) for the two- 
dimensional XY model. Tc/J = 1.05 (red solid line), 1 (green 
dashed line), 0.95 (blue dash-dotted line), 0.9 (dotted pur- 
ple line) and 0.85 (black solid line). The arrows indicate the 
merging points with the line of quasi-fixed points. The verti- 
cal line shows the value of pg. 



and 



Vk 



.92 



92 



-2 


"i + ^c7[ 




L 92 J 



1 

8^ 



(61) 



where we evaluate the rhs in (|6ip at the minimum (which 

we denote by po,/c for simplicity) of the potential Uk{p) 
[Eq. (|37|) ] (see Appendix lA)) . /i, I2 and 1^ are defined in 
([53| . We use the change of variables (|48| with n — i and 

tc = -2. 

The flow trajectories in the plane {po,r]) are shown in 
Fig. [T21 The flow diagram is reminiscent of the results 
obtained in the continuum 2D 0(2) model [s^ [s^l- At 
low temperature {T < J), the trajectories join a line of 
quasi- fixed points. The value of po at the merging point 
depends on the temperature. The critical temperature 
of the BKT transition is defined by the trajectory for 
which the merging point corresponds to pg. A precise 
determination of the value of Pq (which can be obtained 
by fitting the beta function kdkPo.k to j;{po^k ~ Po)^^^) 
is however not possible in the LPA' as it requires the 
full 0{d^) expansion of the effective action in the con- 
tinuum limit fc ^ 1 1131 • Since for the long-distance 
properties of the XY model the lattice should not mat- 
ter, we make the assumption that the ratio PQ/pmax 
(Pmax is the value of po for which i] is maximum) takes 
the same value in the XY and continuum 0(2) models. 
We can then deduce the value of pg from the results of 
Ref. [13. Although this determination of Tc is clearly ap- 
proximate, we can nevertheless conclude from our results 
that 0.9 < Tc/J < 1. This estimate should be com- 
pared to the exact result Tf^^-ct _ q gg j obtained from 
Monte-Carlo calculations [38| and the mean-field expres- 
sion 

T^MF ^ 2J. Note that the relative error on Tc is of 
the same order of magnitude as in the 2D Ising model 
(SecHIEI. 



V. CONCLUSION 

We have proposed a new implementation of the NPRG, 
which takes as a reference system the local limit of de- 
coupled sites rather than a system where fluctuations are 
frozen. The lattice NPRG captures both local and crit- 
ical fluctuations in a non-trivial way. For a lattice fleld 
theory and classical spin models, the LPA is sufficient 
to compute non-universal quantities (transition temper- 
ature and magnetization) to an accuracy of the order of 
1 percent. We have also discussed an approximation (the 
LPA') which goes beyond the LPA and allows to compute 
the anomalous dimension rj. 

A new NPRG scheme has been recently propose d by 
Blaizot, Mendez-Galain and Wschebor (BMW) (H, 
The BMW approach relies on approximate flow equa- 
tions for the effective potential Uk and the two-point 

(2) 

vertex P^ ' . Contrary to the LPA and the LPA', it keeps 
the full momentum dependence of the two-point vertex. 
We believe that the BMW scheme provides the natural 
framework to go beyond the LPA in the lattice NPRG. 

There are many theoretical methods where the idea to 
use a reference system which includes short-range fluc- 
tuations is central. These methods are usually based on 
"cluster" approaches where one solves exactly (usually nu- 
merically) the model on a single cluster (possibly a sin- 
gle site) [4J and includes the coupling between clusters 
by means of a perturbative calculation, a self-consistent 
condition, etc. The cluster approaches include the cor- 
related cluster mean-field theory in classical spin models 
(see e.g. Ref. [l^ and references therein), the dynami- 
cal mean-field theory lisl (DMFT) and its extensions 
(cellular DMFT [isl, |4y|), the cluster perturbation the- 
ory [47I |48| , the variational cluster approach [4^ and the 
self-energy functional theory [sd']. These approaches de- 
scribe exactly the local fiuctuations but struggle to take 
into account low-energy (collective) fluctuations which 
become very important near a phase transition or in low 
dimensions. 

In this context the lattice NPRG may be seen as a step 
towards a theory including both local and critical fluctua- 
tions in strongly-correlated systems. While we have only 
discussed classical models, the lattice NPRG can be eas- 
ily applied to interacting boson systems (Bose-Hubbard 
model) |5i]. 
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Appendix A: Anomalous dimension 

In this Appendix, we compute the anomalous dimen- 
sion for the Ising model. From the deflnition (|45p and the 
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flow equation ^ , one obtains the foflowing expression of 
the (running) anomalous dimension 77^. = — A;9/jhiZj,, 



dg. Gk (q)5,, [kdkRk (q)Gfe (q)^] , (Al) 



where 73 is defined after P5)) . To proceed further, we use 



R = Zekyr, 
kdkR = -{vr + 2yr')Zeky, 
d,^R = Z{r + yr')d,^e, 
dgJdkR = -Z[7jkr + (77 + 4)yr' + 2y\"]d,^e, 



(A2) 



dg^G = -G^Z{l + r + yr')dy^e, 



where r = r{y) = 61(1 - y)(l - y)/y, e = eo(q), G = G(q) 
and y — e/ek- To aUeviate the notations we now drop 
the k index. The product of dq^G oc (l + r + yr') = 9{y~ 
1) and dq^kdkR gives zero except for possible singular 
contributions at ?/ = 1 coming from r"r and r"r' . Thus 
equation (jAip simplifies into 



y\"{l+r + yr')G\dqA' 



(A3) 



By an integration by part we obtain 
y\"{l + r)G\dq^ef 



ek / y\l + r)G^dq^edqy 



where G = Zek/{Zek + f/'(po) + 2pot/"(po)). The second 
contribution in (IA4I) reads 



0(1 -y) 



id,.^)L=,. (A6) 



since 9q^e = 2eo sin((72;) vanishes for g^; tt. For /i — > 
we obtain 



2G^ 



q_L 



where the Lorentzian /i/(q^+/i^) acts as a delta function 
^ '5('?±)- Thus the integral vanishes for d > 2 and takes 
the value G^/(27rZ'*e^.) for d = 2. By a similar reasoning, 
we find 



y\"r'G\dq^e)^ 



G^ 



2{ZtkY 



^(1-2/) 



e(Ce)+ 3(9,^6)2 



(A8) 



We deduce 

7|G4 



77 : 



2t^Z^tl 



k VJ 



i^k - e) 



"d,2 

2tt 



0(qi-M^)/9,^ [y\l + r)G\dq^eY 



'qj- 

.2 ,,2\„,2/ 



^ .ql-A^^)y^(l + r)Gna,,e)r'|;::;, (A4) 



where q^ = (QvtQzi ■ ■ ■ )• Note that we have regularized 
the integrals near q = (/i — > 0+). To compute the first 
integral in (|A4|1 . /i, we remark that if dq_^ acts on G* 
then the integrand vanishes, so that 



(A9) 

In the continuum limit, one recovers the known results 
of the LPA'. In particular, when U{p) = ■|(p — is 
truncated to second order in p, 



77 = 72 



d (l + 2poA)4' 



(AlO) 



G^ 



(ZekY y 



0(1 -y) 



'.<dle) + {dq^ef] 



(A5) 



Note that the last term in (|A9p ensures that 7; is a con- 
tinuous function of d. 

For the XY model, a similar calculation leads to (pDl . 
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